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ABSTRACT 

We have conducted a sensitive search down to the hydrogen burning limit 
for unextincted stars over ~200 square degrees around Lambda Orionis and 20 
square degrees around Sigma Orionis using the methodology of Koenig & Lci- 


sawitz (2014). From WISE and 2MASS data we identify 544 and 418 candidate 


YSOs in the vicinity of Lambda and Sigma respectively. Based on our followup 
spectroscopy for some candidates and the existing literature for others, we found 
that ~80% of the K14-selected candidates are probable or likely members of 
the Orion star forming region. The yield from the photometric selection criteria 
shows that WISE sources with Kg — w3 > 1.5 mag and Kg between 10-12 mag 
are most likely to show spectroscopic signs of youth, while WISE sources with 
Kg — w3 > 4 mag and Kg > 12 were often AGNs when followed up spectro¬ 
scopically. The population of candidate YSOs traces known areas of active star 
formation, with a few new ‘hot spots’ of activity near Lynds 1588 and 1589 and a 
more dispersed population of YSOs in the northern half of the H II region bubble 
around a and e Ori. A minimal spanning tree analysis of the two regions to 
identify stellar groupings finds that roughly two-thirds of the YSO candidates 
in each region belong to groups of 5 or more members. The population of stars 
selected by WISE outside the MST groupings also contains spectroscopically 
verified YSOs, with a local stellar density as low as 0.5 stars per square degree. 
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formation — stars: pre-main-sequence 


HII regions — infrared: stars — stars: 


department of Astronomy, Yale University, New Haven, CT 06511, USA 

department of Astronomy, California Institute of Technology, Pasadena, CA 91125, USA 

3 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 







2 


1. Introduction 

While the data return from the Spitzer Space Telescope vastly improved our picture of 
star formation and circumstellar disk evolution through study of a large number of molecular 
clouds at mid-infrared wavelengths, Spitzer was a pointed mission with a small field of view. 
The more recent WISE mission mapped the entire sky in 3.4, 4.6, 12, and 22 /irn filters to 
a specified 5cr point source sensitivity of 0.08, 0.11, 1, and 6 mJy (or 16.5, 15.5, 11.3, and 
7.9 mag) at spatial resolution of 6-12". With sky coverage better than IRAS and spatial 
resolution within a factor 2-3 as good as Spitzer , a clear advantage of WISE is its coverage 
of many regions of interest that Spitzer simply missed. Sensitivity with WISE is comparable 
at all but the 22 / 24 /irn bands to Spitzer’s shallow surveys, e.g. GLIMPSE in the Galactic 
plane. 

Circumstellar disks surround young stars for the first several Myr of their lives. These 
disks create infrared through millimeter wavelength excesses due to thermal emission from 
dust that is distributed over distances from a few hundredths to several hundred AU from 
the star. Spitzer and WISE both probe the infrared excess emission that arises within a few 
AU. 


As Spitzer mapped primarily the young constituents of known molecular clouds, our 
knowledge of disk evolution time scales based on Spitzer studies is still biased and somewhat 
limited. WISE, however allows us to fill in many gaps in studies of nearby star forming 
regions. Particularly noteworthy wide-held studies that were enabled by WISE include those 
by Rebull et al. (2011) on Taurus, Rizzuto et al. (2012) and Luhman & Mamajek (2012) 
on Upper Sco, and Koenig et al. (2012) and Koenig & Leisawitz (2014) on outer Galaxy 
regions. 


The Orion molecular cloud and young star complex is a benchmark star forming region 
(see Bally 2010 for a review). Orion’s prominence is due to its proximity (~400 pc) as well 
as its representation of both high and low mass star formation, of clustered and isolated 
star formation, and of future, ongoing, and recent star formation. The a Ori cluster, just 
west of the Horsehead Nebula and below the easternmost “belt star” in the constellation 
Orion, is a 3-5 Myr old cluster associated with the central 09.5V star. Its stellar population 


has been studied by Wolk (1996), Sherry et al. (2008), Caballero et al. (2008), and Lodieu 


et al. (2009) and currently contains ~350 confirmed members and another ~300 photometric 


candidates. There is relatively little molecular gas in the region (e.g. Lang et al. 2000); it 
is located between the Orion A and Orion B giant molecular clouds. The A Ori cluster, 
near the famed Betelguese in Orion, is a 4-8 Myr old cluster associated with the central 08 
III star and contains ten additional B stars. Its low mass population has been studied by 


Dolan & Mathieu (1999, 2001, 2002), Barrado y Navascues et al. (2004), and Sacco et al. 
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(2008) though membership is vastly incomplete and knowledge of spectral types limited. A 


surrounding ring includes dust and swept-up neutral and molecular gas (e.g. Lang et al. 


2000) and potentially younger stars than those immediately surrounding A Ori. There is 


also suggestion of a previous supernova explosion (Cunha & Smith 1996). For reviews of 
these two regions, see Walter et al. (2008) and Briceno (2008) respectively, for the a Ori 
cluster and the extended Orion OBla/OBlb region included in the presentation below, and 


Mathicu (2008) for the A Ori region. 


Rich clusters in the age range of cr Ori and A Ori are relatively rare. Complete surveys 
with WISE of these two clusters are thus important for census building, and for overall and 
uniform studies of circumstcllar disk properties and for our understanding of disk evolution. 
Analysis of 2MASS data alone shows that these two clusters indeed have a lower fraction of 
near-infrared excess relative to many 1-3 Myr old clusters still assocatiated with molecular 
gas. A key question is whether the mid-infrared excess fraction is also proportionally lower 
or similar; the former would imply disk depletion at all radii relatively quickly whereas the 
latter would suggest inside-out depletion. 


2. The WISE view of a and A Orionis 


Figure [l] shows a 3-color composite of WISE image data in the 3.4, 12, and 22 /mi 
(wl, w3, w4) bands for the a and A Orionis fields. Emission at 12 /irn (green in the figure) 
is dominated by bright PAH emission that traces cloud surfaces and reveals the well known 
A Orionis ring, as well as the ring around the extended a Orionis field. The 3.4 /irn band 
(assigned to the blue channel in the figure) has some contribution from PAH emission lines, 
but mainly picks out stellar photospheric emission from foreground and background stars. 
The 22 /irn band (in red) captures thermal emission from dust grains which highlights cloud 
surfaces in the same way as w 3, but also heated dust close to OB stars. A SIMBAD search 
shows that in the A Ori region this feature is seen around f 1 Orionis (B0.5I1I) and HR 1763 
(B1V). In the o Ori region a bow-wave feature around a Orionis, (see Ochsendorf et al. 2014) 
and halos around VV Ori (B1V) and HR 1861 (B1IV) are picked out by bright, extended 
22 /irn emission. These stars are labeled in the figure, as well as the three main belt stars of 
the Orion constellation and X Ori, a Mira variable evolved star. 


Note that the other red, evenly spaced, diffuse patches of emission aligned vertically 
with the bright emission in the eastern part of the a Ori region and also appearing in the 
lower, central part of the panel are image artifacts in the 22 /irn WISE images. 
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Fig. 1.— WISE 3-color image mosaics for the a (upper panel) and A Orionis (right panel) 
regions. The image assigns the 3.4, 12, and 22 /im channels to blue, green and red in the 
image respectively. The large, diffuse red spots appearing at regular intervals northward 
from NGC 2023/2024 and from 6 h 35 m , —4° are bright source image latents from wA. 
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3. Initial YSO Selection Criteria 


The combination of 1.2, 1.6, and 2.2 /iin photometry from 2MASS with the 3.4, 4.6, 
12, and 22 /irn photometry from WISE is a sufficient lever arm on the stellar (e.g. from 
J — K colors) to circumstellar (e.g. from wl — w2 and wl — w4 colors) regime that robust 
identification of evolving and dissipating young circumstellar disks can be made. 

We initially extracted sources from the regions surrounding a and A Orionis from the 
WISE All Sky point source catalog. For the former we extracted a region bounded by 
84.0<RA<85.5 degrees and —3.0 <Dec< 2.05 degrees (J2000.0), for the latter a polygon 
with vertices at (RA, Dec) = (77.1 11.6; 87.7 15.4; 90.8 6.5; 80.5 2.8) to encompass the full 
extent of the large ring of emission that surrounds the A Orionis Hll region, as seen in the 
WISE band 3 atlas images at 12 /irn. 

The WISE All Sky catalog also provides photometry in the near-infrared J, H and K$ 
bands from 2MASS. The WISE All Sky photometry pipeline uses a match radius of 3" to 
identify the 2MASS counterparts of WISE sources. 


We then applied an early version of the WISE-2MASS YSO classification and galaxy 
filtration scheme of Koenig et al. (2012). As described by Koenig et al. the scheme first ap¬ 
plies WISE color and magnitude cuts to remove background unresolved star-forming galaxies 
and active galactic nuclei (AGNs). These objects are typically redder and fainter than the 
majority of young stars in Galactic star-forming regions. A further round of cuts is made 
to remove objects likely to be shocked blobs of gas in young star outflows (bright at WISE 
band 2) and spurious detections of bright nebular emission (objects with very red w2 — w3 
color). The remaining sample is then classified into young star ‘classes’ using WISE colors, 
matching the categories laid out by Greene et ah (1994) as best as possible. This scheme 
does not attempt to find either Class III or Flat SED objects, however, only Class II and 
Class I. Additional young stars are selected based on combined WISE and 2MASS colors, 
aiming to retrieve those objects missed due to their WISE band 3 detection being obscured 
by bright nebular background emission. The final steps of the scheme look for candidate 
transition disk objects (Strom et ah 1989) with red w2 — w4 colors and test whether the 
previously classified Class I objects have red enough w2 — w4 colors, or are more likely to 
be the more evolved Class II. 


In the initial selection process for this study, we differed from the Koenig et ah (2012) 
scheme just described in the follow ways. We did not set an upper limit to the magnitude 
error in a given band, but simply required that if photometry was needed for a given clas¬ 
sification criterion, that the photometric uncertainty should have a non-null value. We also 
used a more lenient w2 magnitude cut for star-forming/PAH galaxies of w2 < 13.5 and 
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we did not apply the very last cut of the scheme that removes objects at the blue end of 
the WISE T Tauri star locus. The region of WISE color-color space roughly defined by 
0.5 < w2 — w3 < 1.1 and 0.1 < wl — w2 < 0.4 is occupied by older, more evolved Classical 


Be stars (as subsequently shown by 

Koenig & Leisawitz 

2014 

hereafter, K14), but it can 

also contain young ‘transition disk’ objects (see 

Cieza et al. 

2012 

). Contamination of WISE 


band 4 at 22 /mi by fake detections due to nebulosity produced a large number of candidate 
transition disks in this initial YSO candidate list. We thus required a signal to noise at band 
4 of at least 4.5, and set an upper limit to the w2 — wA of 6 to suppress falsely classified 
objects. 

Having produced a list of candidate YSOs for both a Ori and A Ori, we additionally 
constrained the source lists before each spectroscopic run at the telescope, as described 
below. 


4. Spectroscopic Target Selection and Observations 

In total, 230 sources were observed spectroscopically at the Palomar Observatory 200” 
telescope. All data were collected with the Double Spectrograph (originally comissioned by 
Oke & Gunn 1982 but with many upgrades since that time) using a dichroic at 4800 A and 
a 1200 1/mm red and a 300 1/mm blue grating. Exposures of FeAr and HeNeAr lamps were 
obtained for the blue and red side wavelength calibration. 

The spectra were collected over several years, with evolving selection criteria. 


4.1. 2009 Spectra 


As part of an unrelated program studying low mass stars in cr Ori (Cody & Hillenbrand 


2014), spectroscopic observations were taken in 2009 of known members that were later 


revealed to be of interest from the vantage of the current WISE study. These objects were 
selected for spectroscopy based on their presence in the fields photometrically monitored by 


Cody & Hillenbrand (2010) and lack of a published spectral type at that time. 


As discussed by Cody (2012), observations took place on the nights of 18-21 January 


and 19-20 December with 68 unique low mass stellar and brown dwarf objects in a Ori 
observed. 
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4.2. 2012 Spectra 

For the first run dedicated to WISE follow-up in A Ori and a Ori, we required all 
stars to have non-null photometric error in WISE bands 1 and 4 and to not suffer from 
the diffraction spike, scattered light halo, optical ghost or latent contaminant flags (SQL 
query: cc-flags not matches ’[DHOP]’). To minimize obvious contamination by background 
galaxies, we required the 2MASS Extended Source Catalog proximity column xscprox to be 
either null or > 30". Finally, using the WISE band 4 image atlas we inspected all the YSO 
candidates and rejected objects that were extended or whose 22 /mi emission appeared to 
be offset from the shorter wavelength centroids by more than 2". In both fields we noticed 
a concentration of objects with red colors that appears at magnitudes fainter than J — 16 
and Ks = 14. The spatial distribution of these detections is roughly uniform; thus they are 
likely mostly extragalactic, so we also cut them from further consideration. We obtained 
optical photometry for all sources where available from the NOMAE0 online database and 
restricted our targets to objects with R magnitude less than 18.3. For cr Ori we were able 
to use imaging data from the Sloan Digital Sky Survey]^] to confirm that objects fainter than 
J = 16 are indeed dominated by galaxies. 

In cr Ori we obtained 29 objects with DBSP on 01 January, 2012, UT, selecting those 
YSO candidates with K s < 14 and without existing spectral types in the literature or in 
Cody’s thesis work discussed above. Most objects having J < 11 and wl — w2 > 0.1 were 
observed, as were most objects at all brightnesses with wl — w2 > 0.7. Among objects with 
J > 11, the observations focused on those without previous designations in SIMBAD, i.e. 
WISE is the first that attention has been called to them. A few objects meeting the above 
criteria had roughly equal brightness, close companions within a few arcsec in 2MASS or on 
the telescope guide camera when slewing to the field. These were excluded from observation 
due to concerns about contamination in the excess selection, since we didn’t have enough 
information to attribute w3 or wA excesses between close equal sources. 

In A Ori we obtained 40 objects with DBSP on 02 January, 2012, UT, selecting YSO 
candidates with J < 16 and wl < 14 and a Ks < 14 cut during the observing. We also 
limited ourselves to targets south of Declination +09 so as to avoid the areas previously 
studied by others, namely, the 1 deg 1 2 around A Ori itself and the B30 and B35 regions, with 
the goal being that our work with WISE is unique in the lower Declination areas. 

The prioritization for observations was a combination of brightness (brighter objects at 


1 http://www.usno.navy.mil/USNO/astrometry/optical-IR-prod/nomad 

2 www.sdss.org 



high airmass) and w 1 — w 2 color. Nothing with wl — w2 < 0.15 was observed (i.e., none 
of the Class Ill/transition disk candidates). We focused on the reddest wl — w2 colors. At 
the faint end, some spectra turned out to be active galactic nuclei (AGNs), but in general 
we attempted to observe the faintest YSO candidates. Spatially, the clustering of excess 
objects coincident with the A Ori ring towards the southeast was prioritized; this spatial 
cut also captured many of the reddest wl — w2 objects that were not so faint as to be 
likely extragalactic contaminants. In examining regions to the west of A Ori, almost every 
WISE-selected YSO candidate with J < 12.5 was observed. 


4.3. 2013 Spectra 

In 2013, all targets were in the A Ori field, again with a focus on the western and 
southern parts of the region. We selected YSO candidates with 8 < wl < 13 to avoid both 
infrared-bright asymptotic giant branch stars (AGBs) and faint AGNs, although we did not 
impose a color criterion. In total, 80 sources were observed on 2, 4, 5 February 2013, UT. 

For this run, we did not apply the additional WISE band 4 quality assessment, or the 
requirement on the cc_flags parameter that were used for the 2012 run. One result potentially 
attributable to this modification is that the fraction of observed sources with emission lines 
was lower than in the earlier spectra. 


5. Spectroscopic Data Reduction and Analysis 

Two-dimensional CCD images from both the blue and red channels of the spectrograph 
were corrected for detector bias and flat field and one-dimensional spectra were extracted, 
wavelength calibrated, and flux calibrated, all within the twodspec and onedspec packages 
of IRAF. In 2012, Feige 34, Feige 110, G1 38-31, HD 19445, and HD 199178 were used as 
flux standards while in 2013, Feige 34, HD 93521, Hiltner 600, and G 191B2B served the 
purpose. For our setup, the full spectral range on the blue side was 3525-5075A, and on 
the red side 4660-10,940A, though mediated by the dichroic and the atmosphere. 

The stars were classified visually by comparison to the Jacoby, Hunter & Christian 
(1984) and the Allen & Strom (1995) standards. The DoubleSpec blue side and red side 
spectra were classified independently, with the final spectral type taken as the intersection 
of the ranges concluded on each side. 

Emission and absorption lines were studied in each object. Equivalent widths were 
measured for Ha, Li I 6707 A, the Ca II triplet at 8498, 8542 and 8662 A, and related H 
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and K lines at 3933 and 3968 A. Other emission and absorption lines noted in some spectra 
included [O I] 6300, 6363 A, [S II] 6717, 6732 A, O I 8446 A, Na I 5890 A, He I 5876 A, and 
the Balmer series of Hydrogen out to H5. Several objects with the Balmer jump clearly in 
emission were also identified. 

We found that objects with Kg — w3 > 1.5 mag and Kg between 10-12 mag were most 
likely to show spectroscopic signs of youth such as strong Ha or the Ca II triplet in emission. 
Extremely red objects with Kg — w3 > 4 mag and Kg > 12 were more likely to be AGNs. 


6. Updated YSO Selection Criteria 

The initial WISE-b&sed selection criteria and spectroscopic follow up observations de¬ 
scribed above served to obtain a sample of young stars with optical spectra in the A and 
cr Orionis regions. Subsequent work by K14 showed that this original YSO finding scheme 
generates a large number of spurious candidate young stars. We thus implemented for final 
presentation here the more robust YSO finding procedure detailed in K14. The K14 scheme 
begins in the same way as K12 by removing extragalactic contaminants (AGNs and star¬ 
forming galaxies). It does not include a step to remove detection of sources that appear red 
due to nebular emission because it relies on the fact that YSO candidates occupy loci in 
color space that naturally separate them from these detections. As with K12, the scheme 
then uses WISE and then PITS'A+2MASS color criteria to identify and classify Class I and 
II YSOs. Final steps select candidate transition disk objects with red w3 — w4 colors and 
retrieve candidate protostars from objects flagged as AGNs that are bright at w4. A new 
feature of the K14 scheme is the addition of an explicit phase of AGB star removal that uses 
cuts in the WISE wl — w2 vs. w3 — w4 color-color and w 1 vs. wl — w2 color-magnitude 
diagrams. Table [l] presents a comparison of the number of YSO candidates produced by the 
K12 method described in Section 3 with the number produced by the K14 scheme, over the 
area of sky shown in Figure [TJ The newer selection algorithm clearly finds a greatly reduced 
number of candidates as a consequence of its more conservative design. 

The spatial distribution of YSO candidates found using the K14 scheme in A and cr 
Orionis is shown in Figure [2j To assess the new contribution made by our survey to the 
distribution of YSO candidates in the two regions, we carried out a census of the relevant 
literature. Both the A and a Orionis regions have been surveyed by many authors. The 


main studie 

cr Ori we cc 

s in A Ori are those by 

Hernandez et al. 

(2010 

and Dolan & Mathieu 

2002). In 

>mpilc the Mayrit catalog of Caballero 

2008) 

and the studies by 

Pena Ramirez 

et al. 

(2012 

), Briceno et al. (2005) 

Megeath et al. 

(2012) 

, Bejar et al. (2011 

) anc 

Luhman 

et al. 

(2008 

)• 
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Figure [3] shows how the compiled list of previously known likely members of the two 
regions relate to the newly identified WISE YSO candidates found in this paper. Figure [4] 
shows the distribution of the targets followed up with optical spectroscopy in the two regions. 
The upper panels of Figures [2j [3] and [4] show that our extracted sample of YSO candidates 
extends beyond the boundaries of the A Ori ring. G192.16—3.82 is an HII region at a distance 
of 1.52 kpc (Shiozaki et al. 2011[ ), with an associated young cluster (Carpenter et al. 1993). 
LDN 1617 is known for its association with Herbig Haro outflow sources and is at a distance 
~400 pc (Kajdic et al. 2012). 


6.1. YSO Spectral Slopes 

While the scheme of K14 provides an assessment of the distribution and a simple break¬ 
down of the evolutionary state of the YSO candidates that it selects, a more commonly 
used means of dividing young objects with infrared excess emission is to use the slope of 
their spectral energy distributions (SEDs) in the near and mid-infrared. We follow the 


methodology of Greene et al. (1994) developed for the availability of K-band and 25 /irn 


IRAS measurements, and compute a, the slope of the SED, for each YSO candidate over 
the longest possible wavelength baseline available from 2MASS and WISE (note: we do not 
perform a fit to the SEDs). Ideally we compute a from the slope between the 2MASS Kg 
band to 22 /irn at WISE band 4, requiring non-null Ks band photometric error, signal to 
noise in wA > 5 and 0.3 < y 2 4 < 1-7. If a source does not meet the Ks band requirement we 
derive a from the slope between wl and wA, requiring signal to noise in w 1 > 3. If the band 
4 photometry fails the above criteria we use WISE band 3 with the same test on signal to 
noise and y 2 at w3, paired with either 2MASS Kg or wl as appropriate. Finally we resort 
to the slope between 2MASS H and WISE w2 bands if Kg, w3 and wA all fail the above 
tests, a is computed from the following equation: 


a = 


d log A F x 
dlog A 


( 1 ) 


Note that for all objects, we attempt to remove the effects of dust extinction by looking 


up the value of Ak in the nearest pixel in the A'-band extinction map of Lombardi et al. 


(2011). We compute the extinction in the 2MASS bands using the extinction law of Rieke 


& Lebofsky (1985) and in the WISE bands using the interpolation presented in Koenig & 


Leisawitz (2014). Having derived a value of a for each YSO candidate, we divide them into 


the categories of Greene et al. (1994) as Class I (a > 0.3), Flat (0.3 > a > —0.3), Class II 
(—0.3 > a > —1.6) and Class III (cc < —1.6). Note that the Class III category does not 
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Fig. 2.— Distribution of YSO candidates in a and A Orionis (upper and lower panels 
respectively) as identified by the scheme of Koenig & Leisawitz (2014). Red points and 
triangles: candidate Class I’s from WISE or R / 7S'F , +2MASS. Yellow points and squares: 
candidate Class II’s from WISE or WISE+2MASS. Blue points: candidate transition disk 
objects. 
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Fig. 3.— Distribution of literature YSOs from the papers cited in the text (cyan box points) 
and WISE YSO candidates (red-orange dot points). A, e and a Orionis are plotted as white 
diamonds. The locations of other notable dark clouds or clusters are also marked. 
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Fig. 4.— Distribution of spectra acquired in this paper (green box points), the Cody et al. 
(2010) sample (blue box points) and two YSOs from Riaz et al. (2015) (light triangle points). 
YSO candidates are shown as red orange dot points. A Ori, and in the lower panel e, S and 
cr Orionis are plotted as white diamonds and ( Ori as a black diamond. 
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signify ‘no infrared excess,’ but rather ‘weak excess.’ We tabulate the resulting breakdown 
of source types in Table [lj In both fields, some fraction of the K14 Class II sources are 
categorized as Class III in the a prescription, and some are placed in the ‘Flat-SED’ class, 
although no more than 20% in either case. However, from the K14 Class I sources, 63% 
of the candidates in a Ori move to the Flat SED class under the a calculation, compared 
with 25% in A Ori. In general this suggests that the Class I sources in A Ori extend to much 
redder colors such that their slopes a are consistent with a Class I source in that prescription. 
However, it should be noted that the most embedded clusters in the cr Ori region, those near 
to NGC 2023 and 2024, are poorly characterized in this survey, owing to the high degree of 
saturation in the WISE band 3 and 4 images in that area. This issue limits the robustness 
of the comparison of the a SED class distributions between A and a Ori. 

The Class III category in the a prescription is an overlapping set with K14’s Class II 
and transition disk sources, because its only definition is a steeply negative SED slope. Since 
transition disk candidates are selected on the basis of weak excess in wl — w 2 color and strong 
excess in w3 — w4, most of the transition disk objects in our samples are actually counted 
as Class II YSOs by the a prescription (14/15 from the A Ori sample, 11/13 from a Ori). 


7. Assessment of the WISE-Selected YSO Candidates 
7.1. Spectroscopic Success Rate in YSO Confirmation 

In A Orionis, 75 of our spectra were YSO candidates from the K14 selection scheme. 
In a Orionis, we obtained spectra of 48 K14 YSO candidates. We examined how many of 
these targets showed emission in either the H a (equivalent width < —10 A) or Ca II 8542 
(equivalent width < 0 A) lines as a simple test of ‘youth,’ although it should be noted that 
many M stars show Ha: in emission. In A Orionis, 50/75 objects (67%) met these criteria, 
while in a Orionis 28/48 spectra (58%) met these requirements. While a more general 
indicator of youth is provided by the Lithium Li I 6707 A line, our Palomar spectroscopic 
data are not only too low in resolution for this to be a reliable youth test, but were also 
found to suffer significant uncertainty due to a feature at this wavelength in the flat field 
that can be traced to the long decay time of a dome light, resulting in additional absorption 
at this wavelength in the final spectra. For this paper, we choose not to use the measured 
Li I equivalent width measurements. 
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7.2. Spectral and Mid-Infrared Properties 


In Figure [6] we show color-color diagrams of the K14 WISE YSO candidates in the two 
Orion regions. We overlay the colors of the spectroscopic sample acquired in each held, 
separating those with and without Ha or Ca II 8542 A above our threshold values, and 
distinguish those spectra that correspond to K14 YSO candidates and those that are not 
picked out by that selection process. In the panel for a Ori, we also show the location of 
the two objects investigated by Riaz et al. (2015), both of which are candidate YSOs in our 
selection scheme and both of which appear to be good spectroscopic YSO candidates. 


In both A Ori and cr Ori a significant fraction (28/75 and 21/48 respectively) of spectro¬ 
scopically observed K14 WISE YSOs do not show significant Ha or Ca II 8542 A in emission. 
However, above J — Ks = 1-3, a majority of WISE selected YSOs do exhibit emission in 
these lines. Figure [6] shows that in WISE color-space, a similar increase in the prevalence of 
Ha and/or Ca II 8542 in emission is seen when wl — w2 > 0.3 and w2 — w 3 > 2.2. 


Objects followed up with spectra also exhibit differing distributions of Ha equivalent 
width between the K14 YSO candidates and those spectra not identified by these criteria 


i.e.. 


acquired under either the K12 criteria for infrared excess or as part of the Cody & 


Hillenbrand (2014) sample). Larger EW(Ha) in emission (i.e. more negative) are seen in the 
K14 selected YSO candidates. In Figure [5] we plot the equivalent width in the Ha line versus 
observed WISE wl — w2 color for the full spectroscopic sample, whether matched with the 
K14 YSO catalog (black points) or not (red crosses). In both fields a K-S test showed a 
probability < 10 -4 that the equivalent width distributions were drawn from the same parent 
distribution. The anomalous red cross points with wl — w2 >0.6 were not selected as YSOs 
by the K14 scheme either because they failed to meet the requisite photometric quality 
criteria, or were rejected as AGB candidates because they have very bright wl apparent 
magnitude. 


We conclude that the K14 YSO selection scheme has been tuned relative to the K12 
YSO selection scheme to be an efficient method of identifying Class II and Class I objects 
having significant excess in 2MASS or WISE Liters. The K14 selected YSOs are likely to 
have greater amounts of circumstellar material and therefore perhaps higher accretion rates 
than those objects with only weak infrared excess. 


7.3. Spectra of Objects Outside the K14 WISE Selection 

Approximately half of the optical spectra are of objects not designated as YSO candi¬ 
dates by the K14 scheme. Specifically, in A Ori, 55 of the 130 and in a Ori, 52 of 100 objects 
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Fig. 5. — Equivalent width in the Her line versus wl — w 2 color for the spectroscopic sample 
in the two fields, flagged as a YSO candidate by the K14 scheme (black points) or not (red 
points). 
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Fig. 6.— Infrared color-magnitude diagrams of the K14-selected YSO candidates in A and er 
Orionis (black points, or green points if a spectrum was acquired). Additional plot symbols 
show spectroscopic targets, colored as in figure legend. 
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for which we acquired spectra do not correspond to a K14 WISE YSO candidate. As de¬ 
scribed above, this feature is partly because some of our spectroscopic targets were selected 
with earlier versions of our YSO selection criteria and partly because some came from the 
separately selected targets of Cody et al. The K14 scheme aggessively selects against objects 
likely to have fake detections in the WISE catalog, in particular in WISE band 4, which 
hinders its ability to find objects that have infrared excess only in the longest wavelength 
WISE band. 

As seen in Figure [5j at color w 1 — w2 < 0.2 there is a dearth of infrared-selected YSO 
candidates. Yet some of these objects are moderate Ha emitters. At color wl — w2 > 0.2 
there are a large number of infrared-selected YSO candidates, but still a few objects not 
selected that are moderate Ha emitters, and conversely some objects that are selected by 
K14 but do not exhibit strong Ha emission (those with small equivalent widths). 

In Table [2] we present a summary of the properties of objects not selected under K14 
but for which we have spectra. We list first the total number of spectra that do not match 
to a candidate WISE YSO, then the number that exhibit Ha or Ca II 8542 in emission 
above our threshold equivalent width values (N (young)). We then give the number that 
have red infrared colors, either w2 — w3 > 1, w2 — u4 > 2 or excess color H — Ks > 
— 1.76 x (wl — w2) + 0.9 (N(excess)). Some of these objects not selected under the K14 
scheme have weak or no excess in WISE bands 1 and 2: Wl — W2 < 0.25 which makes them 
candidate transition disk objects, that is, stars with cleared out or depleted inner regions. 
The number of these is given in Table [2] as N(TD — like). Of the full sample of spectra that 
do not match K14 WISE YSO candidates, 24/107 (22%) show spectroscopic signs of youth 
in Ha or Ca II 8542. 


7.4. Contaminants: Spatial Distribution and Colors 

Spectra marked in Figure [6] that appear to be those of AGNs or quasars in the A Ori 
field are a result of the overlap in WISE color space of these extragalactic sources with 
protostars and the reddest Class II sources. These objects escape elimination because of 
their brightness in wl. 

In the A Orionis field, of the K14 YSO candidates with spectral follow-up, four were 
found to resemble those of extragalactic objects, whether AGNs or other galaxies. One K14 
Class I candidate object was found to be a known quasar in a SIMBAD cross-match. One 
object appeared to have a featureless but blue spectrum. In Figure [7] we show a color- 
magnitude diagram summarizing the results of the overlap between the WISE photometric 


Table 1. 


WISE YSO Candidate Breakdown 


Method 

N (Class I) 

N(Flat) 

N(Class 

II) 

N(Class III) 

N(Transition Disk) 

K12 

405 


1621 

A Ori 


12653 

K14 

77 


452 



15 

a 

58 

80 

353 


53 


K12 

159 


854 

a Ori 


2855 

K14 

24 


381 



13 

a 

9 

27 

300 


82 



Note. — K12 refers to the scheme of Koenig et al. (2012), specifically as described in 
Section 3. K14 refers to the scheme of Koenig & Leisawitz (2014). a refers to YSO classes 
based on the slope of the infrared SED. 


Table 2. Spectroscopic targets not selected by WISE scheme 


Field 

Total 

N (Young) 

N (excess) 

N (TD-like) 

A Ori 

55 

10 

35 

22 

a Ori 

52 

15 

6 

0 


Note. — TD-like means ‘transition-disk’ like, as discussed 
in the text. 
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YSO selection scheme and our spectral follow-up program. In the panel for a Ori, we also 


show the location of the two YSOs studied by Riaz et al. (2015) 


The extragalactic spectra were found amongst the protostar candidates in the A Orionis 
held and in general were found in regions of low YSO surface density. In Figure [8] we show 
the distribution of protostar candidate sources color-coded by their projected angular space 
density cr 6 , divided into quartiles of apparent space density. We compute <r 6 by finding 
the distance to each object’s 6th nearest neighbor in the complete YSO catalog (Class II, 


transition disk, as well as protostars) and following the prescription of Casertano & Hut 
(1985). The quartiles are as follows—the lowest 25% of protostars have <Te <0.75 stars 


per square degree, the 2nd quartile 0.75< a 6 <1.89 stars per square degree, 3rd quartile 
1.89< <Tq <5.73 stars per square degree and the 4th quartile have cr 6 >5.73 stars per square 
degree, color-coded from dark to light red. In the same figure we overlay the location of the 
objects whose spectra (or SIMBAD types) were extragalactic or otherwise non-stellar. 


The colors of low density YSO candidates also show they may be largely dominated 
by extragalactic contaminants. In Figure [9] we show color-magnitude diagrams of all WISE 
YSO candidates above and below apparent density 06=0.75 stars per square degree. We 
also show for comparison a contour representation of the color-magnitude distribution of 
stars within ~ 10° of the North and South Galactic poles (requiring w2 signal to noise > 
5, xi < 1-8 an d K s uncertainty < 0.1). The plot of YSO candidates at low spatial density 
shows a relative increase in the quantity of sources with Kg — w3 >4 and Kg >12.5, similar 
to the region in color-magnitude space seen in the Galactic pole reference held data. This 
location is also where the objects found in the Galactic pole reference fields are found. This 
trend fits the expectation that the majority of stars are born in clusters, while extragalactic 
contaminants should be found in a more uniform distribution. 


Koenig & Leisawitz (2014) estimated a contamination rate by extragalactic sources by 


randomly, repeatedly placing an appropriate number density of A11WISE catalog objects, 
drawn from a patch of sky away from the Galactic Plane, behind the 2D extinction map of 


Schlegel et al. (1998) for their Outer Galaxy test held and finding the average number of 


these that would be selected as YSO candidates by their scheme. Using the numbers from 
their Table 12 and scaling for the area of our Orion helds, we can estimate that 41 of the 
77 Class I candidates and 32 of the 452 Class II candidates in A Ori may be extragalactic 
contaminants. In a Ori the numbers are 4 of 24 and 3 of 381 respectively. K14 estimated 


the Galactic contamination rate using the Galactic stellar populations model of Wainscoat 


et al. (1992) to predict an upper limit to the number density of non-YSOs that would appear 


in their test held in the same infrared magnitude and color ranges as YSOs in their scheme. 
We estimate that 28 of the Class I candidates and 93 of the Class II candidates in A Ori 
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Fig. 7.— Infrared color-magnitude diagrams of the K14-selected YSO candidates in a and A 
Orionis (black points, or green points if a spectrum was acquired). Additional plot symbols 
show spectroscopic targets, colored as in figure legend. 
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Fig. 8.— Distribution of protostar candidates in the extended A Orionis field color-coded by 
their local apparent space density. Red boxes: AGNs or quasars, orange box: spectra blue 
and featureless. 
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may be Galactic contaminants and in cr Ori, 3 Class I and 9 Class II’s respectively may be 
Galactic contaminants. Onr spectroscopic sample in A Ori confirms none of the 6 Class I 
sources observed and 40 of the 64 Class II’s as young stars on the basis of Hct or Ca II 8542 
emission line strength. Our spectroscopic sample in a Ori (together with the samples of 
Cody et al. and Riaz et al.) confirms 3 of the 4 Class I sources observed and 26 of the 46 
Class II’s as young stars in the same way. The Class I success rates rate and K14 predicted 
contamination rates appear to be consistent for the two fields. For the Class II sources, we 
confirm only the smaller subset of our candidates that possess Ha or Ca II 8542 in emission, 
and miss out on those objects that may have Lithium in absorption. A more precise test of 
the contamination rate predictions would require a measurement of the Lithium absorption 
line strength and a larger spectroscopic sample. 


The analysis of Koenig & Leisawitz (2014) showed that the YSO candidate lists gener¬ 
ated by their selection scheme were likely contaminated at a rate of roughly 173 extragalactic 
sources and at most 290 Galactic sources in their ~480 square degree test field in the Outer 
Galaxy, or 0.36 sources deg~ 2 and 0.6 deg -2 respectively. These estimates predict 72 extra- 
galactic and 121 Galactic sources contaminating the A Orionis WISE YSO list and 7 and 
12 sources in the a Orionis list. The number of extragalactic contaminants was estimated 
in K14 by randomly, repeatedly placing an appropriate number density of A11WISE catalog 
objects, drawn from a patch of sky away from the Galactic Plane, behind the 2D extinction 
map of Schlegel et al. (1998) and finding the average number of these that would be selected 
as YSO candidates by the K14 scheme. The Galactic contamination rate was estimated by 
using the Galactic stellar populations model of Wainscoat et al. (fl992) to predict an upper 
limit to number density of non-YSOs that would appear in that part of the sky in the same 
infrared magnitude and color ranges as YSOs in the K14 scheme. Using the spectroscopic 
data we have for the two fields combined (including the samples of Cody et al. and Riaz 
et ah), for those sources selected by K14 as Class II YSOs, we End 66 out of 110 objects 
followed up with spectra show either Ha or Ca II 8542 in emission. For the Class I sources, 
3 out 10 objects show either Ha or Ca II 8542 in emission and we hnd 6 likely extragalactic 
contaminants, whether AGNs, quasars or featureless spectra. The total predicted contam¬ 
ination rate using the estimation of K14 above is ~22% for the two fields combined. This 
rate appears low for the Class I sources, albeit from a small sample. The success rate of 
confirmed Class II candidates in these fields appears low as well; however, the absence of 
lithium absorption equivalent width measurements for these objects means there may be 
more young stars within this sample than we can ascertain from emission lines alone. 


We note that Assef et al. (2013) quote a true AGN (dusty quasar) surface density of 
62 per square degree, based on the simple selection criteria w 1 — w2 > 0.8 and w2 < 15.05 
in high latitude fields. This density of objects is much higher than the general density of 
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all YSO candidates in the two Orion fields in this paper. In a Ori the median local space 
density is 33 deg -2 , while in A Ori it is 9.5 deg -2 . The AGN surface density of Assef et al. 


(2013) is also higher than the predictions of Koenig & Leisawitz (2014) for the extragalactic 
contamination rate. K14 estimated that their YSO candidate list was likely contaminated 
by approximately 173 extragalactic sources in their ~480 square degree test field in the 
Outer Galaxy, or 0.36 sources deg -2 . The high surface density of AGNs found by Assef 
et al. (2013) is largely due to their deeper survey in w2. Our YSO candidate sample is 


entirely brighter than w 2 = 13.5, resulting from a cut in w 1 introduced in the K14 selection 
scheme specifically to try to minimize the number of AGNs included in the YSO candidate 
sample. The extragalactic objects in our spectroscopic sample are simply the bright tail of 
the galaxy population. In general, a w 1, w2 or K magnitude upper limit is a critical criterion 
for avoiding extragalactic objects in a spectroscopic survey. 


7.5. Ancillary Information on Candidates 

Tables [3] and [4] present the entire WISE-2MASS derived catalogs of YSO candidates for 
the cr and A Orionis fields respectively. We give the A11WISE catalog designation, coordi¬ 
nates, 2MASS and WISE photometry, source YSO class determined by the K14 scheme and 
the computed SED slope a. We list the equivalent width measurements in the Ho and Ca II 
8542 A lines and spectral types where available from our spectra. 

As a final step in vetting the list of mid-infrared selected YSO candidates, we consulted 
Palomar DSS and SDSS images to determine whether sources could be readily identified as 
galaxies, and also SIMBAD^] for previous information in the literature regarding possible 
status as a contaminant or a young star. Based on this information we have designated 
sources in Tables [3] and [4] with a ”y” (yes) or ”m” (maybe) or ’hr” (no) to indicate the 
likelihood of membership in the Orion star forming region. Where spectral types were 
available in the literature we have included these in the tables and noted the specific reference 
in the final column. 

In A Ori, of 545 K14-selected candidates, 349 have no additional membership informa¬ 
tion from either the literature or our spectra. Of the 196 with supplemental information 
beyond just the 2 MASS + WISE photometry, we designate 149 as probable members (in¬ 
cluding 50 from our spectra) and 26 (including 17 from our spectra) as possible members 
for a confirmation rate of 76-89%. In a Ori, of 418 K14-selected candidates, 147 have no 
additional membership information from either the literature or our spectra. Of the 271 with 


3 http: / / simbad.u-strasbg.fr 
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Fig. 9.— Color-magnitude diagrams of all A Ori WISE YSO candidates above (left panel) and 
below (right panel) apparent space density ct 6 =0.75 stars per square degree. The background 
contours shows a 2D histogram of objects drawn from an equivalent area of sky around the 
North and South Galactic poles in the A11WISE catalog (200 square degrees). Contour levels 
log-spaced starting at 1, increasing in intervals of 0.7dex. Colored points are as in Fig. [7] 
upper panel. 
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supplemental information, 233 are probable members (including 38 from our spectra) and 
33 (including 6 from our spectra) are possible members, resulting in a confirmation rate of 
86-98%. Only 20 sources in A Ori and 4 in cr Ori are ruled out as galaxies or evolved stars. 


8. Final YSO Census 


8.1. Completeness 


A summary of the number of YSO candidates found using the K14 scheme is presented 
in Table [TJ Assuming a 400 pc distance, the K s < 14 limit of the selected YSO candidates 
illustrated in Figure [ 7 ] corresponds at an age of 1 Myr to a mass of 0.05 M 0 , with the mass 
limit increasing as the stars get older (about 0.15 M 0 at 10 Myr and 0.5 M 0 at 1 Gyr). In 
the presence of even modest extinction, we should be sensitive to the vast majority of young 
star candidates having infrared excess. 


To assess the completeness of the K14-derived YSO sample, we take as a comparison 
sample the updated list of young stars found in nearby star forming regions from the Cores 
to Disks (c2d) survey produced by Hsieh & Lai (2013). We take their SpiAer photometric 
catalog, classify their young stars with the scheme of Gutermuth et al. (2009) and produce 
lists of Class I, II and transition disk young stars. We then search for that subset of objects 
in the A11WISE catalog and classify them following K14. At an apparent magnitude of 
K s — 13.25 (10.3) the retrieval rate of Spitzer Class I or II YSOs by WISE is down to 50% 
(90%). An apparent magnitude of A%=13.25 converts to a mass limit of 0.1 M 0 (1.4 M 0 ) for 
an age between 3-5 Myr, using the stellar evolution models of Siess et al. (2000), assuming 
a distance of 400 pc and zero extinction. 


8.2. YSO Spatial Distribution 


In Figure [2j we show the distribution on the sky of YSOs found using the scheme of 
Koenig & Leisawitz (2014). Figure [4] shows the distribution of the spectra we acquired for 


this paper. 


We note in Figure [9] that not all of the low density YSOs are contaminants. Several 
YSO candidates in the right panel have Ks < 12.5 and 0.8< Kg — w2 <2, a locus that is 
not a feature of the reference field color-magnitude distribution. These are likely young stars 
that either have drifted away from their birth clusters or that formed in relative isolation. 


We assess the clustering properties of the YSO candidates in Orion by constructing 
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a minimal spanning tree from the source distribution. An MST is defined as the network 
of lines, or branches, that connects a set of points together such that the total length of 
the branches is minimized and there are no closed loops. We follow the methodology of 


Gutermuth et al. (2009) to identify a branch length cut off that separates groups or clusters 


from the low density background. 

We fit two line segments to the cumulative distribution function of MST branch lengths: 
a steep-sloped segment at short spacings and a shallow-sloped segment at long spacings. We 
then adopt the MST branch length of the intersection point between the two lines as the 
critical cutting length. For the purposes of this analysis, we define a cluster as a group of 
five or more points all with branch lengths shorter than the cutoff. We run this analysis 
on both regions, using the complete YSO candidate lists and show the resulting clusters in 
Figure 10 In A Ori, 65% of the YSO candidates belong to MST-identified groups of 5 or 


more members, while in a Ori the fraction is 69%. 


In the right panel, showing the A Ori region, the MST analysis picks out several stellar 
aggregates away from the A Ori ring. To the extreme north-east (RA > 89°, Dec > 16°), 
it identifies a YSO cluster associated with the more distant G192.16 Hll region. Along the 
southern part of the panel, small groups near to L1617 are seen (RA > 87°, Dec < 6°), 
as well as the northern edge of the larger 25 Ori/Orion OBla cluster (81 < RA <83°, Dec 
< 4°). Within the A Orionis ring our catalog of YSO candidates resides in clusters around 
the head of the B35 pillar (85 < RA < 88°, 8<Dec < 10?5) and around B30 (RA ~ 82?7, Dec 
~ 12?5). In between these groups the MST algorithm divides the YSOs into small groups 
around the central A Ori cluster (purple, pink and green groups). The majority of the ring 
of emission seen in WISE band 3 does not appear to possess a large quantity of young stars 
with infrared excess, however small clusters of Class II and I sources are coincident with 
L1588 and L1589 (magenta and red points, 79 < RA < 81?5, 5<Dec < 9°). 

In the right panel showing the a Ori region, the a Ori cluster itself is identified by 
the MST algorithm (red points, centered on RA ~ 84?7, Dec ~ — 2?7), as are the groups 
associated with NGC 2024 and 2023 (orange points, RA> 85?2, —2.6 <Dec ~ —1?5). The 
more distributed population of YSOs to the north and west of these regions (83<RA < 85°, 
—2.6 <Dec < —0?2) is also picked out, but broken up by the MST algorithm into several 
groups. We are uncertain about the nature of the remaining, small groups in the right panel 
of the figure without further spectroscopic surveys to probe to lower masses and find and 
young stars without disk excess emission. 


Using our spectroscopic sample, we confirm 12 YSO candidates outside MST-defined 
groups in A Ori and 1 outside the MST-defined groups in a Ori. The lowest local stellar 
density around such objects, as characterized by <j 6 (see Section 7.4), is 0.5 stars per square 
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Fig. 10.— MST groups (colored points) and all other YSO candidates (black points) in a 
and A Ori, left and right panels respectively. Stars are colored according to their group for 
visual identification. 
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degree in A Ori and 18.4 stars per square degree in a Ori. This diffuse population of YSOs 
could be made up of stars that have dispersed from where they formed. As with the remainder 
of the young stars we identify though, a complete census of the mass function and the non- 
disked population is necessary to determine the exact nature of these objects. 


9. Discussion 

9.1. The Wide-Field distribution of YSOs in a and A Ori 

While WISE lacks the sensitivity and resolution to match the depth of previous YSO 
surveys in Orion, Figures [2j [3] and [4] demonstrate its ability to capture the basic distribution 
of star formation on large areal scales. 


9.1.1. a Ori region 


In and around a Orionis itself, several extensive and deep surveys of YSOs have focused 
on mapping out the low mass end of the initial mass function, resulting in the large number 
of literature sources marked in Fig. [3] In the brightest region of nebular emission in the field, 
the Spitzer survey of Megeath et al. (2012) also identifies a much greater number of YSO 
candidates than our present WISE survey, roughly 370 YSOs deg -2 , compared to 84 YSOs 
deg -2 from WISE in that location, or a factor ~4 greater. In the remaining northern and 


western parts of the region, the spectroscopic study of Briceno et al. (2005) and our YSO 


candidates comprise a more diffuse population of YSOs. As can be seen in Figure [2j the 
WISE YSO candidates are mostly Class II or transition disk objects. The mismatch between 
the literature YSOs and our sample in this part of the field is due to the ability of Briceno 
et al. to find weak-line T Tauri stars (WTTs) that lack significant infrared excess which our 
color-criteria cannot easily pick out. We find only 2 of the 105 objects noted by Briceno et al. 
as WTT objects in this region, but 21 of their 33 ‘CTT’ objects in a 2" cross-match between 
catalogs. As discussed by Briceno et ah, these northern YSOs are a part of the Orion OBlb 
association. We do not note an obvious preferential association of our YSO candidates in 
this area with the cloud boundary/the bright rim traced by w3 emission however. 
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9.1.2. A Ori region 


In A Ori our catalog of YSO candidates traces the previously observed distribution 
of YSOs that extends from the head of the B35 pillar across to B30, passing through the 
central A Ori cluster. Our catalog captures a subset of the YSO distribution mapped out 


by the Spitzer survey of Hernandez et al. (2010) around the central A Ori cluster. We also 


miss a fraction of the YSOs documented by the works of Dolan & Mathicu that lack strong 
infrared excess. The small cluster associated with L1588 and L1589 appears to be new to 
this study. Previous work in this area has focused on the several Herbig-Haro objects and 
outflows (HH114, 115, 328 and 329) concentrated around IRAS 05155+0707 (see for example 
Connelley et al. 2007). Our catalog also finds a low density distribution of candidate Class 


I sources to the east and north-east of A Orionis, but still within the greyscale mosaic area 
in Figure |2} Two of these objects have already been found to be background AGNs (see 
Section [70]). As noted by Koenig & Leisawitz (2014), the WISE Class I candidate selection 
is known to be the most highly contaminated by extragalactic sources. Because of their lack 
of obvious association with stellar groupings or bright nebular emission, we suspect that 
most of this group are likely background contaminants. 


10. Conclusions 


We have conducted a sensitive search down to the hydrogen burning limit for unextincted 
stars using WISE over ~200 square degrees around Lambda Orionis and 20 square degrees 
around Sigma Orionis. We used the methodology of Koenig & Leisawitz (2014) that builds 
on the heritage of the work of Koenig et al. (2012), to identify a sample of 544 stars in 
Lambda and 418 stars in Sigma that are candidate YSOs. 


We conducted optical spectroscopic followup of 14% of these candidates in Lambda and 
11% in Sigma Orionis. On the basis of strong emission in the Ha or Ca II 8542 A lines, we 
confirm 41 and 27 young stars respectively in the two regions. More reliable spectroscopic 
data covering the Li I 6707 A line would enable us to confirm or reject a greater number of 
the WISE selected YSO candidates. Based on our followup spectroscopy for some candidates 
and the existing literature for others, we found that ~80% of the K14-selected candidates 
are probable or likely members of the Orion star forming region. 


We improved our understanding of the yield from the photometric selection criteria in 
various ranges of color-color and color-magnitude space, in particular WISE sources with 
Ks — w3 > 1.5 mag and Ks between 10-12 mag were most likely to show spectroscopic 
signs of youth. WISE sources with Ks — w3 > 4 mag and Ks > 12 were often AGNs when 
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followed up spectroscopically. 

While we improved the census of young stars in the A and a Ori regions, the presence 
of several strong Ca II and/or Ha emitters among spectra obtained of stars not selected 
via the K14 methods suggests that the criteria are—though robust and reliable—not 100% 
complete in finding the young stars that are present. 


The population of newly identified and candidate YSOs roughly traces the known areas 
of active star formation, but we identify a few new ‘hot spots’ of activity near Lynds 1588 
and 1589. We confirm the more dispersed population of YSOs seen in the northern half of 
the Hll region bubble around a and e Ori surveyed in Briceno et al. (2005). External to the 
A Orionis ring we note clusters of YSO candidates around the L1617 cloud and the more 
distant G192.16—3.82 Hll region. 


We present a minimal spanning tree analysis of the two regions to identify stellar group¬ 
ings. We find that roughly two-thirds of the YSO candidates in each region belong to groups 
of 5 or more members. Given the likely rate of contamination of the YSO candidate sample 
by galaxies and older stars, as described in Section [70 we suspect the fraction of YSOs in 
clusters to be higher than this value. The population of stars selected by WISE outside the 
MST groupings does contain spectroscopically verified YSOs however, with a local stellar 
density as low as 0.5 stars per square degree. 
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